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Abstract 

We consider tiie nonleptonic and semileptonic decays of Dg-mesons into (j) and /o(980) mesons. 
QCD sum rules are used to calculate the form factors associated with these decays, and the cor- 
respondig decay rates. On the basis of data on Df vr+Tr+vr^, which goes dominantly via the 
transition Df — > 7r+/o(980), we conclude that there is space for a sizeable light quark component 
on /o(980). 
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I. INTRODUCTION 



The interpretation of the nature of the hghtest scalar mesons have been controversial 
since their first observation over thirty years ago. Due to the complications of the non- 
perturbative strong interactions there is still no general a^eement about their structure. 
Actually, the observed light scalar states are too numerous [1] to be accomodated in a single 
qq multiplet, and therefore, it has been suggested that some of them escape the quark model 
interpretation. It is not known whether there is necessarily a glueball among the light scalar, 
and whether some of the too numerous scalars are multiquark or some meson-meson bound 
states, or even admixtures of quarks and gluons j^. 

In particular, the structure of the__jneson /o(980) has been extensively debated. It has 
been interpreted as an ss state 

hadrons . . 

The recently measured relative weight of the reaction /o(980)7r"'" —>■ tt^tt'tt^ jlol |. 

may serve as a tool for the estimation of the ss component of the meson /o(980). As a 
matter of fact, if /o(980) has a pure strangeness component (/o(980) = ss), the dominant 
— s> /o(980)7r+ decay proceeds via the spectator mechanism, as shown in Fig. 1. However, 
in the four quark scenario (/o(980) = ss{uu + dd) /2) the decay — > /o(980)7r+ is expected 
to proceed through a much more complicated recombination. 



, as an four quark ssqq state [Oi]? bound state of 
, |8| , and as a result of a process known as hadronic dressing 




FIG. 1: Schematic picture of the spectator mechanism for the decay Df /o(980)7r+. 

Since the spectator mechanism provides a strong production of the 0(1020) meson in the 
decay D'^ 4>'^^, in this work we consider the ratio 

_ TjDt ^ /o(980)7r+) 
F(D+^07r+) ' 

to evaluate the importance of the ss component in the /o(980) meson. This same ratio was 



evaluated in recent calculations by using the spectral integration technique 



and the 



constituent quark meson model [12] • In both calculations the authors concluded that the ss 
component dominates the /o(980) meson and, therefore, the spectator mechanism dominates 
the — s> /o(980)7r+ decay. Here we use the QCD sum rules method to evaluate the ratio 
in Eq. (Q), as well as the branching ratios for the nonleptonic (pir^ and semileptonic 

— s> (pi'^i'e and /o(980)£^z/£ decays. The two branching ratios involving the meson 

will be used to check the reliability of the method, since these two branching ratios are 



known experiment aly 



^expp+ ^ ^^+^^^ _ (2.0 ± 0.5)% . (3) 



II. DECAY WIDTHS 



The decay width of the nonleptonic process Mtt^, where M stands for the or 

/o(980) mesons, is given by: 

TiDj M7r+) = -^^\M\'^X{ml^,ml,ml) , (4) 



IGvrm 

J-' s 

with X{x,y, z) = + y'^ + — 2xy — 2xz — 2yz. The QCD factorization formula (in the 
limit ml 0) gives for /o(980): 

\M{Di /o(980)vr+)p = ^\V^s\W (ci + |)' /.^« - ^JV^W, (5) 
where /+ is the Dg /o(980) form factor defined as 

(/o(p')I^7m(1 - l^MDsip)) = i{f+mp + p'), + f-{t)q,), (6) 
with t = and q = p — p'. And for we have 

|X(D+-07r+)|2 = ^|KsnK.p(ci + |)'^^A(m^^,m|,m2)[(m^^+m^)^(0) 

- (mz,,-m^)A2(0)]^ (7) 

where Ai and A2 are Dg ^ (p form factors defined as j3] 



+ ^^^etp q, + (8) 
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In Ens. ® and (171) the coefficients Ci and C2 are the Wilson coefficients entering the 
effective weak Hamiltonian evaluated at the normalization scale n |14l |: 



'ci(/.) + ^lo, + J, (9) 
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with Oi = {u'j^{l — 75)^) (57^(1 —75)0). Therefore, in calculating the ratio in Eq. ([Q) we are 
free from the uncertainties in the Wilson coefficients and in the CKM transition elements. 
In the case of the semileptonic decays — > Mi^ui the differential decay rates are given 

by 

for _D+ — > /o(980)£"'"z/f. The decay rate for the decay D'^ (pi'^i^i is written in terms of the 
helicity amplitudes 

H±{t) = {mn^ + m^)A,{t) T ^ '''' ^' V (t) , (11) 

Ho{t) = ({ml^ -ml- t){mn. + m,)Mt) - ^("^^-'"^j- 'h ,{t)) , (12) 



so that 



dt lQ2Ti^nvU ^ ^0 



t\'l\m'm)A)\H^{t)\\ (13) 



dVi G\ I Ks P 



-tX^I\ml^,ml,t)\H,{t)\\ (14) 



III. SUM RULES 



The meson in the initial state is interpolated by the pseudoscalar current 

jD,(a;) = s(x)275c(x) , (16) 

where c and s are the fields of the charmed and strange quark respectively. Summation over 
spinor and colour indices being understood but not indicated explicitly. The final hadronic 
state M is interpolated by the current 

iM{x) = s{x)Tms{x) , (17) 
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where 

T.J' fo. M ^ Mm (^3) 

[ 7o for M = 

Using the QCD sum rule technique the form factors in Eqs. © and (jH)) can be 
evaluated from the time ordered product of the two interpolating fields in Eqs. (fTBj) and (fT7|) 
and the weak current = (1 — 75)0 

TMp^p') = ^ I d'xd'y (0|T[jM(x)j7(y)jl(0)]|0)e^(^'-«-^) . (19) 



In order to evaluate the phenomenological side we insert intermediate states for Dg and 
M, we use the definitions in Eqs. dH)) and (jH]), and obtain the following relations 

^'^^"m, + ms (p^ - mljip'^ - 771%) 
+ contributions of higher resonances , (20) 

for /o(980), and 

- -^^ cy^: ip^ - ^ijip- + -^)^-(^)^--" 

+ ■ ■ ■ ^ + contributions of higher resonances , (21) 

for 0, where we have shown only the structures important for the evaluation of the form 
factors Ai, and V . 

In the above equations the coupling of the /o(980) to the scalar current js = ss, was 
parametrized in terms of the constant as: 

{0\rss\fo) = mj, ff, , (22) 

and we have used the standard definitions of the couplings of Dg and with the correspond- 
ing currents: 

(0|jz.J/^.) = , (23) 

rric + mg 

(O|s7c,s|0) = /^e„ , (24) 

The three-point function Eq. ()19|) can be evaluated by perturbative QCD if the external 
momenta are in the deep Euclidean region 

p < (m^ + rrisf, p'^ < Ami, t < (^e + m^)^ . (25) 



In order to approach the not-so-deep-Euchdean region and to get more information on the 
nearest physical singularities, nonperturbative power corrections are added to the pertur- 
bative contribution. In practice, only the first few condensates contribute significantly, the 
most important ones being the 3-dimension, (ss), and the 5-dimension, (sggCr-Gs), conden- 
sates. For each invariant structure, i, we can write 



][ POO 



ds 



du ■ 



Pi{s,u,t) 



+ 7;^='(ss) + Tl'=^{sgsa.Gs) + ■■■ . 



(26) 



The perturbative contribution is contained in the double discontinuity pi. 

In order to suppress the condensates of higher dimension and at the same time reduce 
the influence of higher resonances, the series in Eq. (j26p is Borel improved, leading to the 
mapping 



(p2 - m2)" ^ (n- 1)! (M2)" 



(27) 



Furthermore, we make the usual assumption that the contributions of higher resonances are 
well approximated by the perturbative expression 



POO 



47r2 



ds du 



Pi{s,u,t) 



luo ^"(s-p2)(«-p/2) (28) 
with appropriate continuum thresholds sq and uq- By equating the Borel transforms of the 
phenomenological expression for each invariant structure in Eqs. fl20|) . ()21|) and that of the 
"theoretical expression", Eq. ()26p. we obtain the sum rules for the form factors (at the order 



Cf,e-Os/''''e-V^' U{t) 



so 



4:71^ J{mc+ms)^ Jo 



UQ 



ds / du 



+ 



{-ss) 



-ml/M^ 



—rric + 2m, + 



mliTLs 
2M2 



+ {sgsCr.Gs)e 



s/M-e-^/M'-p^{s,u,t) 



ml{mc — rris) 2mc — 



8M4 



6M2 



■ml{Amc — Srris) — 2t{mc — rris) rric — 2ms m^mg — 2t{mc — mg 



24M2M' 



6M' 



24M2M' 



(29) 



C^{mn.+m,ye-<s/^'\-"'V''' Mt) 
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-1 r^o 



X pi(.,t.,t)] + %^e-^/^^^ 



t — m^ 



-m^ms 



47r2 . ' {rric+ms) 

m^ms 



UQ 



ds / du 



■^-s/M^^-u/M' 



2M2 



(t - mi) 



+ {sgs(y.Gs)e 



1 m2(m2 + 2mcms — t) 2m1 + ^imcm^ — 2t 
'6 ^ 8M4 12M'2 



3m2 + dmcmg ~ 4t ^ mc{2ml + 3mcms + tms)/2 — t{2m1 + 3mcms/2 — 2t) 



12M2 



6M2M' 



(30) 
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— 1 /"^o 



ds 



Uo 



' {m.c+ms)'^ Jo 
1 - 



du 
2M2 



{sgsCr-Gs)e 



X 



6M2 8M4 



+ 



2m^ — rrirm^ 



2t 



12M2M' 



(31) 



and 



-1 /-^o 
47r2 J{mc+ms)'2 

+ (ss)e-™'/^^' (^1 
1 



ds 



UQ 



du 



e-'/''\--/''' pvis,u,t) 



{sgsa.Gs)e 



X 



2M2 

2mi — mf.m^ — 2t 



where 



C 



fo 



rric + rris 



3M2 ' 4M4 ' 6M^M'^ 



and 



(32) 



(33) 



The decay constants /d^, //q and defined in Eqs. ((221), (I2S1), and fj^ . and appearing in the 
constants C/^ and C^, can also be determined by sum rules obtained from the appropriate 
two-point functions. The explicit expressions for the two-point sum rules and for the double 
discontinuities in Eqs. (j^ . (plj). (jM]) and (jH^ are given in Appendices A and B respectively. 



IV. EVALUATION OF THE SUM RULES AND RESULTS 

In the complete theory, the form factors on the right hand side of Eqs. fl29p . (jHUj) . (jHTj) 
and (j32|l should not depend on the Borel variables and M'^. However, in a truncated 
treatment there will always be some dependence left. Therefore, one has to work in a region 
where the approximations made are supposedly acceptable and where the result depends only 
moderately on the Borel variables. To decrease the dependence of the results on the Borel 
variables , we take them in the two-point functions at half the value of the corresponding 
variables in the three-point sum rules Q^Ql- We furthermore choose 



— 9 = ■ 34 

M'2 mlj ^ ' 

If the momentum transfer, t, is larger than a critical value tcr, non-Landau singularities have 



to be taken into account 



Since anyhow we have to stay away from the physical region , 



i.e. we must have t -C {mc+m^Y , we limit our calculation to the region —0.5 < t < 0.4 GeV. 
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1.6±0.1GeV in Eq.L^ ^ ^ 

(35) 

2.0 ± 0.1 GeV in Eqs.dnHI), (ED) and 



In this range the t- dependence can be obtained from the sum rules directly. It can be fitted 
by a monopole, and extrapolated to the full kinematical region. 

Since we do not take into account radiative corrections we choose the QCD parameters 
at a fixed renormalisation scale of about 1 GeV^: the strange and charm mass = 
0.14 GeV, rric = 1.3 GeV. We take for the strange quark condensate {ss) = 0.8{qq) with 
{qq) = —(0.24)^ GeV'^, and for the mixed quark-gluon condensate {sgsO-.Gs) = ml{ss) with 
ml = 0.8 GeV^. 

For the continuum thresholds we take the values discussed in the Appendix A: sq = 
7.7 ± 1.1 GeV^ and, 



We evaluate our sum rules in the range 4.5 < < 9.0 GeV^, which is compatible with 
the Borel ranges used for the two-point functions in Appendix A. In Fig. 2 we show the 
different contributions to the form factors /+, Ai, A2 and V at zero momentum transfer, 
from the sum rules in Eqs. ()29|) . (jHUj) . (jM^ and ()32|) . as a function of the Borel variable 
M^, using the continuum thresholds sq = 8.8 GeV^ and uq = 1.6 GeV^ or 2.0 GeV^ for /o 
or respectively. We see that ^2(0) gets a big contribution from the quark condensate, 
while the perturbative contribution is the largest one for all other form factors. Such kind 
of behaviour had been already obtained in the D K* semileptonic decay studied in . 
The mixed condensate contribution is negligible for all four form factors, and the stability 
is quite satisfactory in the Borel range studied. Varying the continuum threshold Sq in the 
range discussed in Appendix A, and also evaluating the sum rules using or the expressions 
given in Appedix A for the meson decay constants, or its numerical values, we get for the 
form factors at t = 0: 

0.40 < /+(0) < 0.48, 0.71 < V{0) < 0.89 , 

0.32 < Ai(0) < 0.42, - 0.43 < ^2(0) < -0.37. (36) 

The value obtained for f+fO) is smaller than the value obtained for the same form factor in 

n 

ref. [12] by using the constituent quark meson model. 

The t dependence of the form factors evaluated at = 7GeV^ is shown in Fig. 3. 
In the range —0.5 < t < 0.4 GeV^ no non-Landau singularities occur for our choices of 
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FIG. 2: Various contributions to the OPE of the form factors as a function of the Borel parameter 
M^. SoUd curve: total contribution; long-dashed: perturbative; dashed: quark condensate; dot- 
dashed mixed condensate contribution. 

the continuum thresholds. The QCD sum rules results can, in this t-range, be very well 
approximated by a monopole expression 

m 



t ' 



(37) 



for all four form factors. The reult of the fit is also shown in Fig. 3. the different values for 
the pole mass, Mp, for the different form factors are given by: 



(1.6 ± 0.2) GeV for /+(t) 
(4.2±0.5)GeV for Ai{t) 
(8.0 ± 2.0) GeV for A2{t) 
(1.95 ±0.10) GeV for V{t) 



(38) 
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FIG. 3: t dependence of the form factors. The sohd, dot-dashed, dotted and dashed hnes give the 
monopole parametrization, in Ea. (|37|) . to the QCD sum rule results for —A2 (circles), Ai (triangles 
up), /+ (squares) and V (triangles dow) respectively. 

In the case of A2 we get a very hight Mp showing a very weak t dependence. This 
approximate t independence stems from a mutual cancelation in the sum rule of an increase 
in the perturbative and a decrease in the quark condensate contributions. Even for Ai the 
t dependence is much weaker than for /+ and V. It is also interesting to notice that Mp"^ 
is of the same order than the one for the semileptonic D — >■ K*ii'i found in jl3l , a nd Mp"^^ 
is compatible with the ones found for the Dg decays into 77 ^| and D —>■ niui isl ]. 

Having the form factors we can evaluate the decay widths for the (f)7i^ and — >■ 

(pi'^h'e decays, given by Eqs. (jlj) and (|THjl respectively. We obtain the following branching 
ratios: 

B{D+ ^ 07r+) = (2.8 ± 0.7)% , (39) 

and 

B{Df = (1.4 ± 0.4)% (40) 

where we have used Vud = = 0.975, = 0.132 and the values Ci(mc) = 1.263 and 
C2(^c) = —0.513, corresponding to the results for the Wilson coefficients obtained at the 
leading order in renormalization group improved perturbation theory at /i ~ 1.3 GeV 
The errors in the above results were estimated only by taking into account the uncertainties 
in the form factors in Eq. pUj) and should be understood as limiting values for the branching 
ratios. 
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We see that, within the errors, our results are compatible with the experimental results 
given in Eqs. ^ and 0- Therefore, in the case of the Dg decay into (p, we can say that 
the factorization approximation works well. From this it seems reasonable to suppose that, 
if /o(980) has a dominant ss component, the factorization approximation should also work 
well for the Dg decay into /o(980). Using Eqs. (jlj), (0) and ((Tj), and the values for the form 
factors at t = given in Eq.(|H^ we get: 

rW -> /o(980).^) ^ 

In the recently measured spectra from the reaction —>■ tt+tt+tt" jl^, the relative 
weight of the channel '7r+/o(980) is evaluated to be 



BjPt ^ 7r+/o(980))ig(/o(980) 

B{D+ TT+TT+TT-) 

and the ratio of yields 



TT+TT 



(56.5 ± 6.4)% , (42) 



^ +^ = 0.245 ± 0.0281°:°}^ , (43) 
is measured. Taking into account the results in Eqs. and (jlSj) one gets: 



^ TjDt ^ /o(980)TT+) ^ 0.140 ±0.046 

r(D+ ^ 0TT+) B{fo - TT+TT-) • ^ ^ 



Using the branching ratio i3(/o(980) tt+tt^) ~ 53%, the authors in ref. |ll| have 
estimated the ratio i? to be i? = 0.275(1 ± 0.25). In a different way E791, using couple- 
channel Breit-Wigner function jl^, found a non significative qk, that means indirectly a 
non significative contribution for the decay channel /o(980) KK. Thus if we assume 
that S(/o(980) ^ tttt) ~ 1 which implies ;B(/o(980) tt+tt") ~ 2/3 (2/3 being the isospin 
factor), using this in Eq. (ji^ we get 

^<^-:»">^ 0.210 ±0.069. (46) 

^{Dj 07T + ) 

Therefore, from our result in Eq. (PT|) . we conclude that there is a significant deviation from 
the factorization approximation for the /o(980)Tr+ decay. This could be an indication 

that there is a sizeable nonstrange component in the /o(980) meson, or even that the /o(980) 
structure is more complex than indicated by the naive quark model. 



It is interesting to notice that the result obtained in ref. 12| for the ratio in Eq. is 



very similar to our result in Eq. (|^T|) . However, the authors in 
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12 1 concluded that their 



mixture of uu, dd. We believe that this conclusion was reached because 
™si.e.p.eted the experimental .esu.t 0. I„ the. wo.ds. the 



result supports a description of /o(980) as a ss state with a possible virtual KK cloud, but 
with no substantial 
the authors in ref. 

collaboration measured R = 0.62, with a very small error. Since from the Particle Data 
Group we only know that i3(/o — > vrvr) is dominant without knowing the exact number, 
there is still an indetermination in the ratio Eq. ()44j). As explained above, if we consider 
B{fo vr+TT^) ~ 2/3, we arrive at the result in Eq. which is smaller than our result 
in Eq. (HH), leading us to an opposite conclusion compared with 

One possible way to test if there really is a sizeable nonstrange component in the /o(980) 
is through the measurement of the semileptonic — > /o(980)£'*"i/ decay, since in this decay 
we do not have problems with the factorization approximation. Our predicition for the 
branching ratio obtained from Eq.ijTU)), by supposing /o(980) as a ss state, is: 

B{D+ fo{980)i+u) = (0.55 ± 0.10)% . (46) 

Any significant deviation from that will definitively imply in a sizeable nonstrange compo- 
nent in the /o(980) meson, which could be or not accomodated in the naive quark model. 
Therefore, we urge the experimentalists to search for this decay. 



V. SUMMARY AND CONCLUSIONS 

We have presented a QCD sum rule study of the Dj" decays to final states containing 
and /o(980) mesons. We have evaluated the t dependence of the form factors f+(t), Ai(t), 
A2{t) and V{t) in the region —0.5 < t < 0.4 GeV^. The t dependence of the form factors 
could be fitted by a monopole form and extrapolated to the full kinematical region. The 
axial- vector form factors Ai and A2 have a much weaker t dependence than the form factors 
/+ and V. 

The form factors were used to evaluate the branching ratios for the decays Df (f)TT~^ 
and (pi^i^e and we have obtained a good agreement with experimental data. Since 

the evaluation of the decay width, in the nonleptonic decay, is based on the factorization 
approximation, our first conclusion is that the factorization approximation works well in the 
case of the decay D'^ ^tt"*". 

We have also evaluated the ratio D'^ /o(980)7r"'"/(D+ (f)7i^) and we got a result 
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bigger than estimate based on experimental data. Based on the fact that factorization 
approximation works well in the case of the decay 07r+, this result can be interpretated 

as an indication that there is a sizeable nonstrange component in the /o(980) meson, or even 
that the /o(980) structure is more complex than indicated by the naive quark model. This 
hypothesis could be tested by the measurement of the semileptonic /o(980)£+z/ decay, 

since there is no problem with the factorization approximation in the semileptonic decays. 



Acknowledgements: We would like to thank H.G. Dosch for giving us the idea for this 
calculation and F.S. Navarra for fruitful discussions. This work has been supported by CNPq 
and FAPESP (Brazil). 



APPENDIX A: TWO-POINT SUM RULES 



In ref. j3] the two-point sum rule for the /o(980) meson was evaluated by considering /o 
as a ss state. They got: 



u 



{ss) 3 + 



M2 



{sgsQ-.Gs) 
M2 



7712 

2M2 



(Al) 



Considering in the interval 1 < < 2 GeV^ uq = 1.6 ± 0.1 GeV^ they got 

ff, = (0.180 ± 0.015) GeV. 



(A2) 



If we consider = 0, the result for //^ does not change significantly, and we get //^ 
(0.19 ±0.02) GeV. 

n 

The sum rule for is given by [19|: 

r2-ml/M^ 



1 r ..^ (^+^^'y-^^' ,-wM^ I '^^s{ss) I (a'G")^ ^ (A3) 



Considering at most linearly and using uq = 2.0 ± 0.1 GeV^ we get 

/0 = (0.232 ± 0.010) GeV, (A4) 

in the interval 1 < M"^ < 2GeV^, in a very good agreement with the experimental value 
/;"^ = 0.234 GeV jl|. 
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For fi)^ the two point sum rule is given by: 



so 



-rric + — + 



TT J (mc+nia)'- 

rris msmf.\ 



mdsgsCr.Gs) 



1_ 



2M2 



2 ■ M2 y 2M2 
Considering at most linearly and using sq = 8.8 GeV^ we get 

fo, = (0.22 ± 0.02) GeV, 

in the interval 2.3 < < 4.5 GeV^, in a good agreement with the value quoted in ref. 
and also with recent lattice determination |2]j]: fo, = (0.252 ± 0.009) GeV. 



(A5) 



(A6) 



APPENDIX B: PERTURBATIVE CONTRIBUTIONS TO THE THREE-POINT 
FUNCTIONS 



In all this work we take into account the mass of the strange quark at most linearly. 
We have checked that the contibution of terms proportional to and higher powers are 
negligible. The perturbative contributions for the sum rules defined in Sec. Ill are: 
3 

~ A3/2(s,M,t) 

— {2ml — s — t + u){su + rricmsis — t + u))| 6{s — sm), 



p+{s,u,t) 



|n 2mcms{2ml — s — t + u) + rnl{s — t + u) + s{—s + t + u 



(Bl) 



pi{s,u,t) 



2A3/2(s,n,t) 



^rricU A(s, u, t) + 2m^(m^ — s — t + u) + 2st 



+ rris 



-t^ + 2ut^ + ut{2ml -u)- 2mlu{ml 



u 



- s2(3t + 2u) + s(3t^ - 2tu + u{2ml + m))] } 9{s - sm) 



(B2) 



P2(s,M,t) 



2A5/2(s,M,t) 

+ m. 



|mcM( 



s — t — u) 



A(s, u, t) + Qst + 6m^(m^ — s — t + u) 

/ + - t^M + 2m^u^(3m^ + u) - ut^{l{)ml + u) 
s^(4t + 3m) + tM(6m^ + 8m^M + ^2) + ^a^gt^ - + M(2m^ + 3m)) 
s{4£^ - 5t\ - Uu{2ml + u) + u{Qml + 4m^M + ^2))] | ^(^ _ g^^)^ (B3) 



Pv(s,M,t) 



2A3/2(s,M,t) 



|mcM (2m^ — s — t + 



u 



s^ + t^ 



ut 



2uml - s{2t + n) } e{s - sm 



(B4) 



14 



where 
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